CONSTRAINTS ON THE EXTRAGALACTIC INFRARED BACKGROUND 
FROM GAMMA-RAY OBSERVATIONS OF MKN 501 



Todor Stanev 

Bartol Research Institute, University of Delaware, Newark, DE 19716 

Alberto Pranceschini 

Dipartimento di Astronomia - Universita' di Padova, Vicolo Osservatorio 5, 1-35122 Padova, Italy 

ABSTRACT 

We use the new results of the HEGRA detector on the TeV 7-ray emission from 
MKN 501 to set upper limits on the energy density of the cosmic infrared background 
(CIRB). Contrary to previous interpretations of the 7-ray spectrum of MKN 421 as 
showing an intergalactic absorption cutoff at 5 TeV, the observed spectrum of MKN 
501 extends beyond 10 TeV and appears to be unattenuated by 77 collisions with the 
low-energy CIRB photons. The upper limits on the CIRB intensity - derived both 
assuming an a priori shape for the CIRB spectrum and without model-dependent 
assumptions ~ are thus quite strong and come almost in conflict with the observational 
evaluations based on deep surveys of extragalactic sources in the near- and mid-IR. If 
spectra at TeV energies for extragalactic gamma-iay sources like this for MKN 501 
will be confirmed with improved statistics, we may be forced to conclude that the 
process of 77 interaction in the intergalactic space is more complex than expected and 
the average intergalactic magnetic field extremely weak (B < 10^^^ G). 

Subject headings: infrared: general — gamma rays: observations — scattering 



Introduction 



Cosmic history from the decoupling (z = 1500) to the epoch of lighting of the first luminous 
sources (at redshifts 2: ~ 3 to 5) is one of the biggest unknowns of present-day observational 
cosmology. High redshifts and dust extinction during early active phases both degrade the 
energetic optical-UV photons emitted by massive stars, decaying particles, or more exotic energy 
sources, to the infrared wavelengths. A fundamental information on the total energy budget 
associated with astrophysical processes occurring at high redshifts is then provided by observations 
of the cosmic background at infrared wavelengths (CIRB). 

Unfortunately, the infrared domain presents various levels of difficulty to the observational 
astronomer, because of the huge backgrounds from the Earth's atmosphere, the Interplanetary 
dust (IPD), and diffuse dust in the Milky Way, in addition to the background produced by 
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the telescope itself. Also the sensitivity and stability of infrared detectors are far poorer than 
those used in the optical. Because of all this, the detection and characterization of the diffuse 
background flux of low-energy photons coming from primeval structures has been exceedingly 
difficult so far. Even dedicated experiments exploiting cooled platforms outside the atmosphere, 
among which the most important is DIRBE on COBE ( [Hauser 1996| ), have failed so far to detect 
significant signals from the CIRB above the intense foregrounds. The upper limits allowed by 
the foreground emission deconvolution of DIRBE maps are significantly higher than expectations 
over a substantial - and crucial - A-range from a few to ~ 100 //m. The situation, at these 
wavelengths in particular, is not likely to improve in the future, until a mission flying to the outer 
Solar System will get rid of the fundamental limitation set by the IPD (both the scattered light 
and dust re-radiation). 

Under these circumstances, a very interesting alternative to the direct detection of CIRB 
photons has been suggested by Stecker and de Jager (1993) soon after the discovery of high-energy 
photon fluxes coming from distant Blazars (with GRO, [Hartman et al. 1992 and with the 



Whipple Observatory, Punch et al. 199^ ). The idea is to infer the CIRB spectral intensity from 



combined GeV and TeV observations of a set of active galactic nuclei (AGN), by exploiting the 
7 — 7 interactions and pair production between the AGN high-energy photons and low-energy 
background photons in the line-of-sight to the source. The interaction is expected to produce an 
absorption feature, testable in principle, in the source TeV spectrum. Interesting limits have been 
discussed by Stecker and de Jager (1993), de Jager, Stecker and Salamon (1994), and Dwek and 
Slavin (1994), all based on TeV observations of the Blazar MKN ^21. 

To summarize, the best current upper limits on the CIRB in the spectral range from 10 to 
40 micron are those reported by de Jager, Stecker & Salamon (1994), with a 2a upper limit of 
A/a < 2 10"*^ W/m^/sr. In the same large waveband interval, marginal detections, at levels of 
2 10-9 < A/a < 2 l{)-^W/m'^/sr were reported by de Jager, Stecker & Salamon (1994) and Dwek 
& Slavin (1994). At shorter wavelengths, 1 to 10 /im, upper limits have been obtained by Stecker 
& de Jager (1993), Stecker (1996), Dwek & Slavin (1994) and Biller et al. (1995). The most 
conservative bounds, accounting for the precise spectral shape of the CIRB, given by Dwek & 
Slavin (1994, A/a < 10^^ W/m? /sr), still keep a substantial factor (>10) higher than the expected 
contribution of known sources. 

While all previous analyses relied uniquely on TeV observations of the blazar MKN 421, we 
exploit here a new high-quality dataset of TeV gamma-ray observations by HEGRA of MKN 501 
during a state of high activity (Aharonian et al 1997a) to further constrain the intensity of the 
diffuse IR background. 

Two sets of constraints on the CIRB spectral intensity are derived in Section 2. One is 
based on the assumption that the background spectrum is dominated by the contribution of 
distant galaxies, and thus reflects the average galactic IR spectrum. The other constraint is free 
of model-dependent assumptions and treats the CIRB as a combination of twelve bins wherein 
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the background spectrum is flat (in A/a) with independent arbitrary normahzations. Extremely 
tight constraints on the CIRB ensue from this analysis, reflecting the missing evidence of any 
absorption in the TeV spectrum of MKN 501. A discussion is given in Section 3 in terms of a 
low average past emissivity both of galaxies and of primeval energy sources. Indeed the limits 
are so severe that they begin to conflict with the integrated IR flux of distant galaxies, recently 
detected in large numbers by ground-based and space observatories. We finally emphasize the 
alternative possibility that the interaction of high-energy gamma-rays with low-energy photons is 
more complex than previously supposed. Hq = 75 km/s/Mpc is assumed throughout the paper. 



2. Data analysis and results 

We use the data set taken in March and April 1997 with the HEGRA stereoscopic system of 
four imaging Cherenkov telescopes ( Hermann 1995D . At that time MKN 501 was in an extremely 



high state, the far brightest gamma-ray source in the sky. The flare was observed by all TeV 7-ray 
observatories in the Northern hemisphere ( Breslin et al. 1997| ). The high gamma-ray flux allowed 



the HEGRA observatory to measure the gamma-ray energy spectrum in short time intervals. The 
spectrum measured between March 15 to 20 covers 1.2 orders of magnitude in photon energy, 
from 0.8 to 12.6 TeV, divided into twelve logarithmically spaced energy bins. The photon number 
spectrum at TeV 7-ray energies is consistent with a power law of differential spectral index a = 
2.49ib0.11. This spectrum is shown in Figure 1 (data points) together with the allowed range of 
best-fitting power-laws derived by the observational team (shaded region). Similar power-law 
spectra for MKN 501 have been observed by various other groups (see e.g. Protheroe et al. 1997). 

The main difficulty with the derivation of the CIRB spectral intensity from TeV 7-ray 
absorption is the lack of knowledge of the production spectrum, which may deviate from a 
pure power-law and/or show absorption at the source. Spectra measured within short time 
intervals, during which the source's activity state is not likely to vary, are very valuable in 
this respect, although there is not an evidence for strong spectral variability of the MKN501 
emission( Aharonian et al. 1997b"D . For this reason we confined our analysis to a spectrum based on 



data collected between March 15 to 20, rather than the one using the whole March-April database 
in the updated version of Aharonian et al (1997a). 

As a first step in understanding the observed spectrum, we attempted to derive the allowed 
range of differential spectral indices a at the source and the CIRB absorption by fitting the 
observed 7-ray spectrum with a variety of assumptions for a and for the intensity Im of the 
extragalactic IR background. We first assumed that the source spectrum in the observed range is 
a pure power-law and that the low-energy background is dominated by the integrated contribution 
of distant galaxies. In such a case, the shape of the CIRB is well constrained, and should reflect the 
average galactic spectrum, which has a minimum at A ~ 10/im corresponding to the intersection 
of the stellar with the dust emission component. We refer here to the detailed spectral shape 
estimated by Franceschini et al (1991). 
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Then the fit was performed in a two parameter space consisting of the source gamma-ray 
spectral index a and the CIRB intensity Im normalized to the model spectrum /o(A) by 
Franceschini et al (1991). For a given spectral index a the absorption due to the CIRB was first 
calculated and the resulting spectrum was then normalized to the detected 7-ray flux above 1 
TeV. This procedure enabled us to avoid introducing a third parameter, the source spectrum 
normalization. Table 1 shows the optical depth assuming Iir/Iq = 1, for all twelve experimental 
energy bins. 

Figure 2 summarizes the results of fitting the observed MKN 501 spectrum. It plots contours 
of the 2D distribution in the parameter plane, including the unphysical region corresponding 
to negative Iir/Iq values, where gamma-rays arc 'created' in collisions with CIRB photons. The 
darker an area is the better the fit is (except for fits with < 1.45, corresponding to the blank 
inner region). The formal best fit (x^ = 1.41 p.d.f. for 10 degrees of freedom, confidence level 
0.16) is in the "unphysical" region: a = 2.52, Iir/Iq = -0.3. The best fit in the "physical" region 
occurs at Iir/Iq) = (no-absorption, xt = 1-42 p.d.f.) and yields a = 2.48, in agreement with the 
result of the experimental group. Flatter TeV 7-ray source spectra allow for higher intergalactic 
absorption: a value for the spectral index of a = 2.16, corresponding to Ijr/Iq = 1.70, is the 
hardest spectrum with xt < l-^O (corresponding to the 95% confidence limit). The best-fit spectra 
corresponding to Iir/Iq=0 and the fiattest one assuming Iir/Iq = 1.7 are also shown in Fig. 1. 

It is important to understand that the fit quality, and correspondingly the limits on the CIRB, 
depend crucially on the manner in which the 'theoretical' power-law fluxes are normalized in the 
fitting procedure. The normalization to the observed 7-ray flux above 1 TeV causes very different 
spectral indices to flt the observations equally well and generates a broad x^ valley in Fig.2. This 
would change significantly if the observed 7-ray energy range were extended on either side. An 
observation at > 300 GeV (the threshold energy of the Whipple telescope) would improve 
the overall normalization and allow us to distinguish better between different spectral indices 'at 
production'. An extension to higher 7-ray energy, hence larger optical depths (see Table 1) would 
make easier the detection of any absorption by the CIRB. 

We finally attempted to obtain model-independent limits on the CIRB, with no a-priori guess 
at the background spectrum. To compute them we conservatively assumed that the high-energy 
source spectrum is the flattest allowed by the flts of Fig. 1, i.e. a = 2.16. We normalize the 
spectrum at the source to be F-y + 1.64 x SFy in the first experimental bin. This normalization 

factor is a source of uncertainty. It is partially justified by the fact that MKN 501 is close enough 
(136 Mpc) compared to the optical depth (600 Mpc for A/a = 7.64 10"^ W m"^ sr~^) that the 
content of this bin would be absorbed only by 20%. 

The absorption of 7-rays of energy E TeV is due to IR photons within a certain wavelength 
interval around the maximum absorption at Cmax = '^{meC^)'^ / E^ {^max = ^■'^'^/i^maxj^V) /im). 
We have then assumed that the absorption of any given 7-ray energy bin is caused by CIRB 
with flat AJ(A) spectrum. The envelope of all upper limits to the CIRB intensity obtained in this 
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way is shown as a histogram in Fig. 3. The hmits are assigned within wavelength intervals that 
contribute 90% of the optical depth for that bin. As we see, the limits become less stringent at 
longer wavelengths. The limit in the 3.4 fim to 24 /im range comes from one single energy bin (5.0 
to 6.3 TeV) where a particularly large 7-ray flux was measured. 

Note that the obtained upper limits fall already very close to recent direct evaluations of the 
IR background based on deep IR surveys ^Franceschini et al. 1991 ). As shown in Fig. 3, the upper 
limits heavily rely on the error in the TeV flux measurement and are less stringent for A > 10 fim, 
where the TeV 7-ray statistics is not as good. 



3. Discussion 

In either case, both assuming the CIRB shape and relaxing it, the constraints on the CIRB 
intensity appear dramatic. Essentially the MKN 501 gamma-ray spectrum does not display the 
expected effect of absorption, it rather shows a simple E~'^'^ power-law spectrum. 

How reliable are these limits in view of the possible systematic errors of ~25% in the energy 
estimates of the HEGRA telescopes ( [Aharonian et al. 1997c )? This is very easy to estimate in the 



case of a flat binned \I\ CIRB spectrum. The limits in Figure 3 would move upward and towards 
shorter wavelengths with the fractional amount of energy overestimate. Similar relaxation would 
occur also in the case of a more specialized model CIRB spectrum. One could use the optical 
depths from Table 1 to estimate the amount of relaxation. Similarly, a higher normalization of the 
source spectrum would relax the model-independent limits by the ratio of the two normalizations. 

Is this featureless spectrum of MKN 501 inconsistent with that observed for MKN 421, which 



is almost at the same distance? The latter has been interpreted by some authors (e.g. Stecker 



19961) as showing a turnover at e ~ 3 — 5 TeV , which was attributed to 77 absorption with the 
CIRB. In fact, new observations of this source during a high activity state do not appear to 
confirm the presence of absorption (Krennrich et al. 1997), and show significant counting rate 
above 5 TeV. So, at least during this high state, MKN 421 seems to show a power law spectrum 
similar to the spectrum discussed here for MKN 501. 

The constraints on the CIRB intensity from TeV observations of MKN 501 start to approach 
the "measured" lower limits at 2.2, 6.7 and 15 /xm given by the integrated emission of galaxies 
already resolved in deep integrations at those wavelengths. Deep surveys have been performed 



from ground in the K-band and from space by the mid-IR camera (ISOCAM, see Cesarsky et al 



19961) on the ISO satellite in the two latter bands. A summary of these "direct" determinations of 
the galaxy contribution to the CIRB, and a discussion of the related uncertainties, are given by 
Franceschini et al (1997) and Oliver et al. (1997). In any case, the CIRB cannot be lower than 
reported at these three wavelengths. 



Few possibilities are left. The first one is that the CIRB is very close to the limits allowed by 
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the gamma-ray spectrum of MKN 501 observed by Aharonian et al. (1997a). This would imply 
a very strong constraint on any signals unrelated to the emission of distant galaxies (see e.g. 



Rowan-Robinson & Carr, 1988, for a review). 



But, in view of the fact that the same power-law spectral shape as show in Fig. 1 for MKN 
501 has been confirmed by later integrations on this source (Aharonian et al. 1997a; Protheroe 
et al. 1997), that apparently a similar shape is also suggested for MKN 421, and because an 
appreciable CIRB flux has already been detected and resolved into discrete sources, we find more 
likely that we have to revise our concepts about the propagation of TeV gamma-rays into the 
intergalactic space, and that something complicates the process. 

The question is: why the propagation of TeV gamma-rays in intergalactic space should not 
produce the expected absorption in high energy spectra of distant sources? A possible solution 
could be that part of the source spectrum is regenerated in 7-ray cascading (pair production -|- 
Inverse Compton). In such a cascading process, the 7-ray spectrum of the source is depleted 
around the region of maximum absorption. If the 7~ray emission of MKN 501 at the source 
extends above 10^^ eV, the spectrum would be depleted in collisions with microwave background 
photons. The e~^e~ pairs generated on the microwave background would Inverse Compton scatter 
on the microwave background to regenerate photons of lower (TeV) energy, thus generating 
'bumps' on a power-law production spectrum ( Protheroe fc Stanev 1993| ). The resulting 7-ray 



spectrum may then appear unattenuated at observation. This would however require not only a 
7-ray spectrum extending to very high energy, but also a very low (~ 10~^^ Gauss) value for the 
extragalactic magnetic field in the direction of MKN501. Otherwise the e~^e~ pairs would deflect 
in the magnetic fleld and form a halo around the source, well outside of the angular resolution of 
the HEGRA detector. 

A good deal of constraints useful to disentangle between these two possibilities are soon 
expected by improved observations of the MKN501 outburst (which has been observed by the 
Whipple and CAT Cherenkov telescopes, Breslin et al. 1997| ) with different energy thresholds and 



wavelength bands and by refined forthcoming data on MKN 421. 
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FIGURE CAPTION 

Fig. 1. — The energy spectrum of the TeV 7-rays from MKN 501 observed by the HEGRA 
detector, together with its power-law spectral fits (lightly shaded area). The solid line shows our 
best fit in the physical region (see text). The dashed line shows the fit with a = 2.16 and Iir/Iq 
= 1.70. 

Fig. 2. — The range of spectral indices and CIRB energy densities that allow for 90% confidence 
fits of the observed TeV 7-ray spectrum of MKN 501, assuming a theoretical shape for CIRB - see 
text. The best formal fit, shown with a cross, is in the unphysical region. The scale is given in 
the upper right corner. 

Fig. 3. — Upper limits for the CIRB density derived from the TeV 7-ray spectrum of MKN 501: 
a) assuming the theoretical shape shown with a thin line; b) assuming flat The thick line 

shows the estimate of Franceschini et al (1991) and the three data points are from Franceschini et 
al (1997). These are three direct evaluations of the CIRB spectral intensity due to faint galaxies at 
A = 2.2, 6.7 and 15 /L^m, based on galaxy models fitting deep counts performed at these wavelengths. 
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Table 1: Optical depths and CIRB wavelength ranges responsible for the absorption of TeV gamma 
rays by MKN 501 for a distance of 136 Mpc {Ho = 75 km/s/Mpc). Column 1 shows the 7-ray 
energy range, columns 2, 3 & 4 give the optical depth, CIRB wavelength of maximum absorption 
and the wavelength range responsible for 90% of the optical depth for the model of Franceschini et 
al (1991). Columns 5, 6 &: 7 give the same quantities for constant XIx = 7.64 10~^ W sr"^. 
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